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^- ■ Abstract. We present here a study based on Xmm-newton data of RX J0256. 5+0006, a medium distant (z=0.36) 

cluster of galaxies found in the Bright SHARC catalog. The X-ray emitting intracluster medium shows a bimodal 
structure: one main cluster component and a substructure in the west, which very likely falls onto the main 
cluster centre. The subcluster shows after subtraction of the main cluster component a cometary shape pointing 
away from the main cluster centre, suggesting that ram pressure stripping is at work. Despite the indication of 
interaction between the two components we surprisingly do not find any sign of temperature gradients, which is 
contradictory to predictions from hydro dynamical simulations of cluster mergers. 

Due to the non-symmetric form of the main cluster we extract three surface brightness profiles in different sec- 
,_^ ■ tors around its centre. We see large variations between the profiles, which we quantify by /3 model fitting. The 

Oh corresponding r c 's vary between 0.1-0.5 Mpc and the /3's between 0.5-1.2. The variations of the profiles and the 

(3 model parameters indicate that the main cluster is not entirely relaxed. This hypothesis is strengthened further 

5_l ■ by the fact that the cluster is over luminous with respect to the (z-evolving) L x — T relation found for nearby 

"^ | clusters. 

Galaxy clusters show a high degree of self-similarity. Comparing our profiles to the scaled reference emission mea- 
sure profile of Arnaud et al. based on nearby clusters, we find that only the profile extracted north-east (NE) of 
the main cluster centre is similar to this reference profile. This indicates that only the NE profile is representative 
S^ ■ for the relaxed part of this cluster component. 

Based on the /3 model parameters of the NE profile and the spectroscopically fitted temperature of kT = 
4-91q4 heV we find for the total mass within rsoo using the hydrostatic approach M500 ~ 4 x 1O 14 M0 for 
the main cluster component. This value is in good agreement with the value (M500 = 3.9 x 1O 14 M0 ) obtained 
using the z-evolving M500 — T relation from the HIFLUGCS sample based on nearby clusters. A non-z-evolving 
M — T relation is only marginally consistent with our result. This is an indication that there exists evolution in 
the M — T relation, as predicted from simple scaling laws. Calculating the corresponding gas mass fraction we 
find f g ~ 18 — 20% which is in good agreement with other work. 

We also develop a simple on-axis merger model for the cluster. As input we use the projected distance of the 
subcluster to the main cluster centre and the velocity difference of main and subcluster based on four galaxy 
redshifts spectroscopically measured with the Kitt Peak telescope. Together with a simple ram pressure model we 
find that the most likely physical distance of the subcluster to the main cluster lies between 0.6 < d < 1.0 Mpc. 
The coherent results of on-axis merger and ram pressure model suggest that the merger in this cluster is indeed 
on-axis and not an off-axis merger with a large impact parameter. 

We find for the ratio of subcluster to main cluster mass values between 20-30% which indicates that the merger 
in RX J0256. 5+0006 is a major merging event. 

Key words. galaxies: clusters: individual: 1. Introduction 

RX J0256. 5+0006 - galaxies: clusters: general - cos- 
mology: observations According to the scenario of hierarchical structure for- 
mation, large-scale structures like galaxy clusters form 

through the accretion of smaller units. Clusters of galax- 

Send offprint requests to: Doris Neumann, ies are still forming today via merger events (Kempner, 

e-mail: dneumann@cea.fr Sarazin & Ricker 2002 ; Neumann et al. 2003). During 
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cluster mergers, the hot X-ray emitting intracluster 
medium (ICM) is heated locally by shock waves and com- 
pression. X-ray observations of the ICM allow to study 
the dynamical state of clusters of galaxies. The spectro- 
imaging capabilities of Xmm-newton and Chandra al- 
low studies with unprecedented precision of the ICM 
in interacting clusters up to relatively high redshift 
(Markevitch & Vikhlinin 2001 ; Worall & Birkinshaw 
2003 ; De Filippis, Schindler & Castillo-Morales 2003). 
These observations contribute to the question on how large 
scale structure and galaxy clusters form and grow, which 
is connected to the evolution of internal physical cluster 
properties as well as cosmology. 

We present in this paper the EPIC Xmm-newton 
observations of RX J0256. 5+0006, which is a medium 
redshift (z=0.36) cluster of galaxies first detected in 
ROSAT PSPC data and is a member of the Bright SHARC 
- Serendipitous High-rcdshift ROSAT Cluster Survey — 
cluster catalog (Burke et al. 1997 ; Collins et al. 1997 ; 
Romer et al. 2000 ; Burke et al. 2003). This cluster is one 
of several SHARC clusters selected for X-ray follow-up by 
Xmm-Newton (see Arnaud et al. 2002b ; Majerowicz, 
Arnaud & Neumann (2002) ; Lumb et al. 2003). 

The paper is organized as follows: in section 2, we de- 
scribe the observation and the data treatment. In section 
3, we present our imaging analysis. Our study shows that 
RX J0256. 5+0006 is most likely composed by two cluster 
components, one main cluster component and a smaller 
subcluster West of the main cluster centre. Sec. 4 focuses 
subsequently on the surface brightness profile of the main 
cluster. This is followed by section 5, in which we present 
a spectro-imaging analysis which comprises a search for 
temperature variations via the construction of a hardness 
ratio map as well as spectral fits. We conduct a spectral 
analysis in distinct cluster regions, selected in function of 
the likely merger scenario found in this particular cluster. 
In section 6 we discuss the main cluster properties. This 
is followed by section 7, in which we present optical obser- 
vations of RX J0256. 5+0006 and compare them to Xmm- 
NEWTON data. Using the previously obtained results, the 
dynamics of the merging event is modeled in section 8. 
This allows us to put constraints on the merger geometry 
and the physical distance between the two cluster compo- 
nents. We also give estimates on the luminosity and mass 
of the subcluster. The discussion and conclusions of all our 
results are given in section 9. 

In this paper, we use a flat and low density cosmol- 
ogy with Ho=50km/s/Mpc, fl m =0.3 and f2A=0.7. In this 
cosmology, at the redshift of the cluster (z=0.36), one arc 
minute corresponds to 423 kpc. Uncertainties and errors 
are 90% level unless specified otherwise. 

2. Data treatment 

2.1. Observation 

RX J0256. 5+0006 was observed with Xmm-newton dur- 
ing its 217-th revolution for a total of 25.3 ks. The 



calibrated event lists were retrieved from the Science 
Operation Centre. We concentrate in this article on EPIC 
data. The thinI filter was used during the observation. 
The observation mode for EMOS and EPN cameras was 
Full Frame and Extended Full Frame Mode, respectively. 

Galaxy clusters are extended X-ray sources with rel- 
atively low surface brightness emission. Thus contamina- 
tion of background emission can play an important role, 
especially in the outskirts of clusters. Therefore proper 
background subtraction is crucial to obtain reliable re- 
sults. As background, we use in the following blank sky 
observations (one for each camera) provided by Lumb 
(2002) 1 . The corresponding files consist of several high 
galactic latitude exposures performed by Xmm-newton 
. To avoid contamination, bright point sources present in 
these observations are cut out. 

For our analysis, we only select single (pattern 0) 
events, for the EPN data and pattern to 12 events 
for EMOS data sets. This selection was performed for the 
source and background observations. 

2.2. Vignetting correction 

The sensitivity of the EPIC cameras is not constant but 
decreases with increasing off-axis angle. To correct for this 
effect, called vignetting, we use the method described in 
Arnaud et al. (2001). In this method each detected event 
obtains an individual weight factor which is the ratio of 
the central effective area at the event energy over the ef- 
fective area at the energy and location of the event in 
the camera. Therefore when creating images, spectra or 
surface brightness profiles, the event weight factors are 
added up into bins instead of the individual event counts 
themselves. Thus, a proper correction of the vignetting is 
assured. The error propagation takes into account the fact 
that we use weight factors and not photons. 



2.3. Background treatment 

2.3.1. Flare rejection 

The background of Xmm-newton shows sometimes pe- 
riods of high intensity which are the result of soft protons. 
These periods are also called flares. To illustrate these 
flares the light curve of the RX J0256. 5+0006 observation 
is displayed in Fig.l. Since these periods of high inten- 
sity limit severely the signal-to-noise of the observation, 
we discard them without further study. 

For the flare rejection, we use the method described 
in Majerowicz, Neumann & Reiprich (2002) i.e. we only 
keep time intervals where there are less than 15 counts 
per 100 s in the 10 to 12keV energy band for each EMOS 
camera and 22 counts per 100 s in the 12 to 14keV energy 
band for EPN (please note we only use Pattern=0 events 



The background event files can be down-loaded from 
ftp://xmm.vilspa.esa.es/pub/ccf/constituents/extras/background/ 
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Fig. 1. Light curve of the RX J0256. 5+0006 observation for the 
EPN camera. Only events with an energy greater than 300 eV 
are selected and subsequently binned in time intervals of 10 s. 



Fig. 2. Combined EMOS1, EMOS2 and EPN image of 
RX J0256.5+0006 in the 0.3 to 2.0 keV energy band. The im- 
age is convolved with a Gaussian (<r=6.6") filter. The black 
contours are the residuals after the subtraction of the elliptical 
13 model from the data (see section 3). The contours are ex- 
pressed in terms of statistical significance (the lowest contour 
is 3 a and the step width is la). The center of the best fit 
elliptical j3 model is represented by the black cross. 



for EPN). The remaining exposure times are 10550s for 
EMOS1, 10340 s for EMOS2 and 7100 s for EPN. 

2.4. Remaining background components 

Beside the particle background there exist two other kinds 
of background: a background induced by highly energetic 
particles which is relatively constant with time and is well 
described by our used blank sky observations, and the cos- 
mic X-ray background (hereafter CXb) which is dependent 
on the sky position (see Snowden et al. 1997). 

To correct for these different background components, 
we use the method described in Majerowicz, Neumann & 
Reiprich (2002) and Pratt, Arnaud & Aghanim (2001). 
We outline it briefly in the following. 

The intensity of the high energy particle background 
varies typically within 5 or 10% from observation to ob- 
servation. In order to account for these variations, we nor- 
malize the background at energies above lOkeV where the 
emission is dominated by these particles. To calculate the 
normalization factor for the EMOS cameras we divide the 
count rate in the 10 to 12 keV band of the source obser- 
vation by the count rate in the same energy band of the 
blank sky observation. Since the EPN camera has a higher 
quantum efficiency at high energies we use for the EPN- 
data the 12 to 14 keV energy band instead of the 10 to 12 
keV band. 

To correct for the local CXB which is not necessarily 
well represented by the blank sky observations, we perform 
a double background subtraction: first, the blank sky ob- 
servations are subtracted from regions inside and outside 
of the cluster. The latter represents the local background. 
For the subtraction we use the normalization factor calcu- 
lated from high energy events mentioned above. Then, the 
remaining residuals of the region outside the cluster are 
subtracted from the blank sky subtracted regions in the 



cluster. The result is thus corrected for all the different 
background components. 

This method will be described with more detail in the 
case of spectra (Sec. 5. 2), surface brightness profiles (Sec. 4) 
and images (Sec.5.1). 

In order to avoid source confusion 46 clearly detected 
point sources are excluded from our analysis. 

3. Cluster Morphology 

Fig. 2 shows the Xmm-newton image of 
RX J0256.5+0006 in the 0.3 to 2.0 keV energy band. The 
cluster exhibits a clear bimodal structure. An important 
question is whether the two components which are 
responsible for the bimodality are physically connected 
and in interaction or whether they are merely due to 
a chance alignment. We will address this issue in the 
following. 

The projected distance <i m i n of the two components 
is about 0.8' which corresponds to a physical distance 
of 350 kpc at the redshift of RX J0256. 5+0006 with our 
adopted cosmological parameters. 

Since the cluster shows an elongation in the North 
South direction we model the hot gas distribution of the 
main component, with an elliptical (3 model (see also 
Neumann & Bohringer 1997). The elliptical (3 model has 
the following form : 



S(x, y) = S (l + F 1 + F 2 )- 3f ' +1/2 + B 



(1) 



with 



Fi = 



F 2 = 



((x — x c ) cos a + (y — y c ) sin a) 
(—(.x — x c ) sin a + (y — y c ) cos a)) 



where (x c , y c ) are the coordinates of the cluster center, M 
is the major axis, m the minor axis, and a the position 
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Parameter 


best fit 


P 


0.83 


M 


370 kpc 


m 


300 kpc 


a 


95° 


■L c 


2 h 56 m 34.3 s 


yc 


00° 06' 11" 



Table 1. Elliptical f3 model best fit parameters of the main 
component of RX J0256. 5+0006. These values are obtained by 
fitting the X-ray emission distribution of the main cluster by 
excluding the western substructure. 



angle. The background level is included via the parameter 
B. 

The fitting technique relics on Gaussian statistics via 
the x 2 test. However, in the outer regions of the image, the 
number of photons per pixel is low and more appropriately 
described by Poisson statistics. To correct for this effect, 
we apply to the image a Gauss filter with a a of 6.6" (see 
also Neumann 1999 for more detail). We also exclude point 
sources in the image and the western structure. The best 
fit parameters are listed in Tab.l. The ratio m/M ~ 0.80 
confirms the apparent cllipticity of the main component. 

Since the emission from the western component is con- 
taminated by the main component emission, we subtract 
the elliptical (3 model best fit (see Tab.l) from the data in 
order to determine the morphology of this substructure. 
The residuals obtained in terms of significance are shown 
in Fig. 2. 

The western structure is very significant at a detection 
limit which largely exceeds 3 a. The residuals can be rel- 
atively well described by an ellipse with major and minor 
axes of 1.7' and 1.4', respectively. At the redshift of the 
cluster, this corresponds to 720 kpc x 590 kpc. These val- 
ues indicate that this structure is very likely a small galaxy 
cluster. Henceforth, we will refer to this western compo- 
nent as the " subcluster " . The morphology of this subclus- 
ter resembles somewhat a comet and indicates strongly 
that this object is in interaction with the main cluster : 
the elongated structure suggests that gas is pushed out of 
the subcluster due to ram pressure stripping as it encoun- 
ters the ICM of the main cluster. 

To examine the robustness of the comet-like shape of 
the structure, we perform two different tests : 

- we vary the (3 value between 0.75 and 0.9 for the ellip- 

tical f3 model ; 

- we subtract a spherical (3 model instead of the elliptical 

model. Two spherical models are used, the first with a 
core radius of m, the second with a core radius of M. 

In each case, the particular comet-like shape of the sub- 
cluster remains. 

The comet-like structure, indicating that ram pressure 
stripping is acting on the subcluster, suggests that the 
subcluster is indeed interacting with the main cluster cen- 



tre. We will discuss the distance between main cluster and 
subcluster in more detail in Sec. 8. 

Numerical simulations (e.g. Ricker & Sarazin 2001) 
show that during a merger event of two clusters, the hot 
gas of the smaller cluster passes through the main cluster 
core once then goes back to the core to be absorbed. In this 
frame, the orientation of the RX J0256. 5+0006 subclus- 
ter isophots indicates that this cluster has not yet passed 
through the main cluster center. Markevitch et al. (2002) 
found for the cluster 1E0657-56 a " bullet " structure com- 
parable to the comet-like structure of RX J0256. 5+0006. 
However, in the case of 1E0657-56, the orientation of the 
isophots revealed that the " bullet " already traversed the 
cluster core. 

Additionally, the residual map also shows significant 
residuals inside the main cluster region. In particular, an 
offset is detected between the maximum of X-ray emis- 
sion and the center of the best fit elliptical (3 model (see 
Fig. 2, black cross and northern excess). This will be more 
discussed in Sec.sbp Sec.7.1. 

4. Surface brightness profiles of the main 
component 

We extract vignetting corrected surface brightness profiles 
of the main cluster. We only select events in the 0.3 to 
3.0 keV energy band to optimize the signal-to-noise ratio. 
We discard regions where point sources are present and 
we apply a binning of 3.3" which is the width of each 
concentric annulus. 

One surface brightness profile is produced for each 
EPIC camera and the so obtained profiles are then subse- 
quently summed in one single profile. We also extract sur- 
face brightness profiles from the blank sky observations 
with the same detector coordinates. These background 
profiles are subtracted from the cluster profiles. The re- 
maining background emission is due to the sky variations 
of the CXB. To correct for this emission, we estimate the 
remaining mean surface brightness value in the annulus 
between 4' and 7' which is well beyond the detected clus- 
ter emission and we subtract this value from the profile. 
We finally group the data into bins with i) a signal-to- 
noise ratio of at least 3c above background and ii) at the 
same time a logarithmic binning in which the binsize at 
radius 9 is: A9 > 1.15 x 9. This kind of binning is similar 
to the binning presented in Arnaud, Aghanim & Neumann 
(2002a). 

In Sec. 3 it was shown that there is an offset between 
the fitted center of the main cluster and the maximum 
of the X-ray emission of roughly 8" . This offset corre- 
sponds to roughly 60 kpc. The location of the maximum of 
the emission coincides with the location of the main clus- 
ter galaxy and might be contaminated by a point source 
such as a central AGN. However, in order to obtain mono- 
tonically decreasing surface brightness profiles, we choose 
the X-ray maximum of the cluster as the profile extrac- 
tion center. We evaluate the potential error based on this 
choice of center by comparing the subsequent Id-fit results 
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Fig. 3. X-ray emission contours from the 0.3 to 2.0 keV image 
of RX J0256. 5+0006 (Fig.2) and the surface brightness extrac- 
tion areas where the blue ellipse region is excluded. 



with the 2d-fit results, which possess different (3 model 
centres: the (3 from the 2d-fit ((3 = 0.83 - please note: we 
do not estimate the uncertainty of this parameter in our 
2d-fit) lies within the found /3-values of the Id-fit of the 
overall profile with and without the extraction of the cen- 
tral excess emission (see Tab.l and Tab. 2). Furthermore 
the geometrical mean of the minor and major axis core 
radius of the 2d-fit (0.3 Mpc), is similar to the core radius 
found for the Id-profile without cutting out the central ex- 
cess region. We conclude therefore that the choice of the 
center for extracting the surface brightness profile is not 
important and that the potentially existing point source 
in the cluster center does not influence significantly our fit 
results. 

Fig. 3 shows the X-ray contours of the cluster. As one 
can see the contours in the North-Eastern (NE) part of the 
cluster are more compressed than in the South-Eastern 
(SE) part. Because of this we extract surface brightness 
profiles in different regions, which are defined in Fig. 3: the 
first region covers the whole area with the exception of an 
ellipse defined in Fig. 3 surrounding the subcluster. The 
second region only represents the NE part of the cluster 
and the third profile is selected in the SE cluster part. The 
resulting profiles are shown in Fig. 4. 

The overall cluster profile can be detected out to a 
radius of 3.5'. This corresponds to a physical radius of 
roughly 1.5 Mpc. Even with our adopted logarithimic bin- 
ning (see above) it is not possible to obtain the chosen 
signal-to-noise ratio (3<r) for radii greater than this value. 

Not surprisingly the NE and SE profiles show very dif- 
ferent shapes. The SE profile is much steeper than the 
NE profile at large radii and the core of the SE profile is 
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Fig. 4. Background corrected surface brightness profiles from 
the main component of RX J0256. 5+0006 in the 0.3 to 3.0 keV 
energy band (the three epic camera, EMOS1, EMOS2 and 
EPN are added). The different regions chosen for the extrac- 
tion are shown in Fig. 3. The triangles show the overall profile 
without angular selection, the circles correspond to the NE 
area and the red boxes to the SE area. The lines show the best 
fit /3 model profiles: full line for the NE-profile, dotted line for 
the overall profile and dashed line for the SE profile. Please 
note: the fits take into account PSF-effects and therefore the 
fitted profile of the NE profile lies above the data points in the 
innermost region. 



much larger than the NE core. This is equivalent to the 
fact that the contours in the NE part of the cluster are 
more compressed than in its SE part. In order to quantify 
the differences, we fit a /3 model to the different surface 
brightness profiles. 

The (3 model gives the relationship of the electron 
density as function of radius (Cavaliere & Fusco-Femiano 
1976) : 



n c (r) = n c0 




(2) 



in which the free parameters are n c o, the central electron 
density, r c , the core radius and [3 a slope parameter. The 
surface brightness S x for the (3 model can be written as : 




-3/3+i 



(3) 



where the S x q is the central brightness. 

We fit the surface brightness profiles with a (3 model 
convolved with the Point Spread Function (hereafter PSF) 
of the detectors (see Ghizzardi 2001; Arnaud et al. 2002b). 
We perform two (3 model fits for each radial profile: the 
first by using all the bins and the second by excluding 
the first five bins to avoid the influence of the detected 
excess, which might be linked to a point source (Sec. 3). 
From the best fit parameters including their errors and 
the mean temperature of the main cluster (see Sec. 5. 2) of 
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Parameters 


Overall 
a b 


Radial profiles 

NE 
a b 


1 
a 


3E 

b 


Xred (dof) 

n e „ (lO^cm -3 ) 
r c (kpc) 


3.3 (22) 

r 9 o+0.20 

286^33 

fl 770 + 6. 059 
u -' '°-0.049 


2.0 (17) 

4.201°^ 

3781™ 

0.8911°:°?? 


1.2 (19) 1.2 (14) 

9.02l°$ 7.02li; 9 i 

1171*°, 153l 8 J 

0.572l°;°f 7 0.6111°;^ 


0.9 (17) 
4.051°;" 

5451^° 

-, ol +0.60 
L - zl -0.26 


0.9 (12) 

4.051$:?* 

5461?^ 

-, 99 +0.74 
1 - zz -0.30 



Table 2. Best fit j3 model parameters for the different extracted surface brightness profiles of the main cluster, a: the fit is 
performed over the whole radial profile, b: the first five bins (corresponding to 120 kpc or 17"- the central excess emission) are 
excluded for the fit (see section3). 



4.9lg'4keV, n c0 can be determined. The results are shown 
in Tab. 2 and illustrated in Fig. 4. 

Not surprisingly, the best fit parameters are different 
for the different profiles. Particularly, one can note that 
the best fit [5 for the SE profile is very high. The fit values 
also change when the central part of the profile is cut out 
by giving generally larger (3 values. However, the differ- 
ent fit parameters found for the entire profile and for the 
same profile with the central region cut out agree generally 
within the error bars. 



5. Spectro-imaging analysis 

5.1. Hardness ratio map 

There are strong indications that RX J0256. 5+0006 is in 
a merger state. In this case it is very likely to observe tem- 
perature variations in the ICM. To visualize these possible 
variations we construct a hardness ratio image of the clus- 
ter, which consists of diving a background subtracted im- 
age in the hard energy band by a background subtracted 
image in the soft photon band. In order to ensure an opti- 
mal choice of the energy bands for the hardness ratio im- 
age we simulate an Xmm-newton spectrum with galactic 
absorption nH=6.6xl0 20 cm~ 2 (Dickey & Lockman 1990), 
z = 0.36 and kT = 4.9 keV (see also below). Based on this 
spectrum we choose for the soft band 0.3-1. 3 keV and for 
the hard band 1.3-7.0 keV. With this choice we ensure to 
have the same number of photons in each band, which 
minimizes the error in the resulting hardness ratio map. 

Since the cluster has low surface brightness and in or- 
der to have sufficient counts in each image pixel to ob- 
tain statistical significant results for the hardness ratio, 
we need to adopt a large pixel size. From the previous fits 
(see Sec. 3 and 4), we compute the total count rate coming 
from the cluster inside a radius of 2.5' in the hard band 
and we find about O.lcts/s. By assuming that the emis- 
sion is flat, we estimate a pixel size of 14.3" x 14.3" in 
order to have at least 3 detected cluster photons in each 
pixel of the hard energy image. 

To correct for the background, we also extract images 
in the same energy bands from the blank sky observa- 
tions. We subtract the background from the cluster im- 
ages according to the normalization defined in Sec. 2. 4. To 
correct for the CXB contribution, we first extract surface 



Fig. 5. Hardness ratio map of RX J0256. 5+0006 obtained by 
dividing the background corrected image in the 1.3-7.0 keV en- 
ergy band by the background corrected image in the 0.3-1.3 keV 
band. Before the division the two images were smoothed by a 
Gauss filter with a=14.3". In the region in which cluster emis- 
sion is found, the temperatures vary between 4.0 and 5.5 keV 
(blue colors correspond to lower temperatures than green col- 
ors.). The contours of the substructure (see Sec. 3) are overlayed 
in black. 



brightness profiles in the soft and hard energy band for 
each camera. Then, we subtract the profiles of the blank 
sky observations from the cluster profiles in the same en- 
ergy band. The residuals at large radii (outside the clus- 
ter emission region) give the correction for the local CXB 
contribution which is assumed to be constant across the 
entire field-of-view. We subtract the obtained averaged lo- 
cal CXB contribution from the blank sky corrected image 
pixels in the soft and hard energy band. This background 
correction is performed for each camera image. The im- 
ages are subsequently summed to obtain a single image in 
each energy band. 

All our cluster images are thus fully background cor- 
rected. To avoid too strong statistical fluctuations from 
pixel to pixel, we finally apply a Gauss filter with cr= 14.3" 
before the division of the hard energy by the soft energy 
image. The final hardness ratio map is presented in Fig. 5. 

Converting the hardness ratio into temperature gives 
values for each pixel which vary between 4.0 keV and 
5.5 keV. This excludes the two pixels in North- West 
which show an extremely high hardness ratio (correspond- 
ing to temperatures greater than lOkeV in the case of 
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Fig. 6. EMOS1 (black), EMOS2 (red) and EPN (green) spec- 
tra of RX J0256. 5+0006 in the region 2 (see also Fig. 7). The 
presented spectra are vignetting and background corrected (see 
sections 2.2 and 2.4). 



brcmsstrahlung) . These pixels correspond to a point-like 
source. 

The hardness ratio map does not show any sign of sig- 
nificant temperature gradients in the cluster. This result 
is surprising for the area between the two interacting clus- 
ters in which heating due to merging is expected. 

5.2. Temperature estimates 

To verify the results given by the hardness ratio map pre- 
sented in Sec.5.1, we perform spectral fits in order to ob- 
tain more reliable temperature estimates. 

Since we apply the described vignetting correction 
method (see Sec. 2. 2), it is possible to use directly the 
on-axis response matrices : ml_thinlv9q20t5r6_all_15.rsp 
(EMOS1), m2.thinlv9q20t5r6_all.15.rsp (EMOS2) and 
epn_ef20_sY9_thin.rsp (EPN) for the spectral fitting. To 
avoid remaining uncertainties at very low energy, we ex- 
clude for our analysis events below 0.3 keV. 

We extract spectra in different regions (see Fig. 7) and 
correct for background (see Sec. 2. 4. We estimate the CXB 
contribution from an annulus between 4.5' and 11'. We 
define the following regions: the main cluster center (Ti), 
the outer parts of the main cluster (T 2 ), the subcluster 
region (T 3 ) and the area in front of the subcluster region 
pointing towards the main cluster centre (T4). 

We group the background subtracted spectra such that 
each bin has a signal-to-noise ratio greater than 3c after 
background subtraction. 

We use Xspec (Arnaud 1996) to fit the resulting spec- 
tra with an absorbed and redshifted isothermal model 
(WABSXMEKAL, see Morrisson & McCammon 1983, Mewe 
et al. 1985, Kaastra 1992, Liedahl et al. 1995) with an ab- 



Region 


T (keV) 


A (solar) 


Xrcd (dof) 


1 
2 

1+2 
3 
4 


4 q+u« 
4 7+ - 8 
4 q+0-5 

r 7 +1.3 

°-'-i.o 

r 9+2.O 
°' z -1.2 


36 +u ' 2y 

U - OD -0.21 
fl 4Q+°33 
u -* o -0.29 

0.36i° ;i? 
0.30 frozen 
0.30 frozen 


1.1 (203) 
0.9 (154) 
1.0 (322) 
0.9 (108) 
0.8 (52) 



- Table 3. Results from the isothermal model fits of the com- 
bined EMOS1, EMOS2 and EPN spectra in the regions defined 
in Fig.7. 



Fig. 7. X-ray emission image from the 0.3 to 2.0 keV image 
of RX J0256. 5+0006 (Fig.2) and temperature estimates in our 
selected cluster regions. 



sorption fixed to the galactic value of nH=6.6xl0 20 cm -2 
(Dickey & Lockman 1990) and a redshift fixed to 0.36. 
The resulting best fit parameters are presented in Tab. 3 
and Fig.7. To illustrate the quality of the spectra, we show 
in Fig. 6 the spectra used for the estimate of T 2 . 

The fact that Ti and T 2 are similar indicates that the 
main cluster does not host a cooling flow. Furthermore, 
we confirm with our temperature estimates (see Sec.5.1) 
the absence of strong temperature variations. The ICM 
appears fairly isothermal. We will discuss this in more 
detail in Sec. 8. 

To check for the proper estimation of the redshift we 
leave this parameter as a free fit parameter for region 1+2. 
The corresponding fitted redshift is 0.351^o'o23 which is in 
agreement with the redshift determined from optical spec- 
tra (see Romcr et al. 2000 and Sec. 7). Since there might 
exist a large velocity difference between the main cluster 
and the subcluster we also perform another fit leaving the 
redshift as a free fit parameter in region 3. Unfortunately, 
the lack of statistics does not allow to constrain the red- 
shift in this region. 

6. Main cluster properties 

6.1. Luminosity estimates 

Due to the pressence of background, we do not detect the 
cluster emission up to the virial radius. It is thus possible 
that we understimate the total count rate of the cluster 
since we miss out on the external parts of the cluster. 
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To overcome this potential problem we use the ellipti- 
cal [3 model fit parameters (presented in Tab.l, based on 
photons in the energy band 0.3-2.0 keV) for the cluster 
emission, which we extrapolate up to the virial radius r v , 
which is defined below. We calculate the corresponding 
count rate of the model extrapolated to r v and convert it 
into a bolo-metric luminosity estimate. 

The case of evolution of galaxy cluster properties with 
redshift is still an open issue since many studies give op- 
posite results (e.g. Vikhlinin et al. 2002 ; Borgani et al. 
2001). Since RX J0256. 5+0006 is not a nearby cluster, 
we choose to apply an evolution term according to the 
self-similar model (Bryan & Norman 1998 ; Eke, Navarro 
& Frenk 1998) to all the scaling laws used in this paper. 
According to the spherical collapse model, the virial radius 



r„ is 



with A z — 



3894A z ^(l + z) 
and \l z = 



T 



lOkeV 

n m (i4 



kpc 



18ir 2 n z 



i l rn \1_ 



a, 



(4) 



where A c is the cluster density contrast (see Bryan 
& Norman 1998) and Vt z the density parameter of 
the Universe (for more detail see Arnaud, Aghanim & 
Neumann 2002a). We find A z = 0.42 for our adopted cos- 
mology and thus r„=2650kpc. This corresponds to 6.3' 
and confirms that we do not detect the main cluster of 
RX J0256. 5+0006 up to the virial radius but only up to 
0.55 xr„ (Sec. 4). However, the missing contribution be- 
yond the detection radius is estimated to be only about 
3% which is similar to the statistical uncertainties. 

For the conversion into luminosity, we use a cluster 
temperature of 4.9lg'4keV (see Sec. 5. 2). We obtain a bolo- 
metric luminosity of Ltoi.x = 1.49 ± 0.10 x 10 45 ergs/s. 

If we apply the measured L x -T relation found by 
Arnaud & Evrard (1999) by assuming a standard evo- 
lution term according to the self-similar model : 



1/2 



L x 



A 



1/2 
z=0 



(1 + ^)210 



45.06±0.03 er S 



6keV/ 



2.88±0.15 



(5) 



the measured temperature of 4.9* ; 4 keV gives a luminos- 
ity of 8. 7+ 2 ; 3 x 10 44 ergs/s (the uncertainties are only based 
on temperature uncertainties, since they are much larger 
than the uncertainties based on the photon number.). 
According to equ. (5), a temperature of about 5.9 keV is 
needed to reproduce a luminosity of 1.5 xlO 45 ergs/s. This 
temperature exceeds the temperature measurements in- 
cluding errors by 0.5 keV. This discrepancy of observed 
and predicted luminosity suggests that this cluster is not 
entirely in equilibrium. 

6.2. The scaled emission measure profiles of the main 
cluster 

Arnaud, Aghanim & Neumann (2002a) found that scaled 
emission measure profiles of distant clusters are similar in 



shape and normalization to the average scaled emission 
measure profile of nearby clusters. 

In the following, we compare the emission measure 
profile of the main cluster of RX J0256. 5+0006 with the 
emission measure profile of Arnaud, Aghanim & Neumann 
(2002a) to sec whether this cluster, dispite its distorted 
morphology follows the found self-similarity of clusters. 

6.2.1. From the surface brightness to the scaled 
emission measure profile 

The surface brightness profile is translated into the emis- 
sion measure (EM) profile via : 



EM(r) oc (1 + z) 



4 S(9) 



(6) 



e(T,z) 

9 is the angle, d\ is the angular diameter distance, 
r = 9 x dj±, and e(T, z) is the emissivity in the considered 
energy band (0.3 to 3.0 keV in our case). The determina- 
tion of e takes into account the response of the cameras 
and galactic absorption(see also Sec. 5. 2). 

The resulting emission measure profile is scaled to the 
self-similar model by using the standard scaling laws and 
the empirical M gas — T relation found by Neumann & 
Arnaud (2001). The physical radius is scaled according to 
the calculated virial radius r v (see equ. 4- please note: this 
value only depends on the measured cluster temperature 
and not on our determined mass profile) and the emission 
measure is scaled by : 



ScEMoc A z s(l + z) ^T 



-1.38 



(7) 



where A z is defined in equ. 4. 

These scaling relations are applied to the three sur- 
face brightness profiles (see Sec. 4). The resulting scaled 
emission measure profiles are shown in Fig. 8. 

6.2.2. Comparison with the local reference profile 

We show in Fig. 8 the scaled emission measure profiles of 
the different cluster regions on top of the reference profile 
of nearby clusters (see Neumann & Arnaud 1999). Only 
the NE profile (see Tab. 4) matches the profile from the 
nearby clusters. The overall as well as the SE profile lie 
systematically too high with respect to the nearby clus- 
ter profile up to a radius of 0.35 x r v . This implies that 
only the NE profile shows self-similarity with other clus- 
ter profiles and at the same time that only the NE profile 
represents the relaxed part of the cluster. The two other 
profiles comprise (a) region(s), which are very likely dis- 
turbed due to past merging activity. The difference of the 
profiles between each other and with respect to the mean 
nearby cluster profile shows how important it is for self 
similar studies to select relaxed (parts of) clusters. 

6.3. Total mass estimates using hydrostatic equilibrium 

In the following we estimate the total mass of the cluster 
using the hydrostatic equilibrium hypothesis. 
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Table 4. Total mass, gas mass and gas mass fraction at rsoo for the main cluster of RX J0256. 5+0006 based on hydrostatic 
equilibrium assuming isothermality. The quantities are calculated using the estimated mean temperature (Sec. 5. 2) and the 
different f3 model best fit parameters (see Sec. 4 and Tab. 2). The uncertainties (90% confidence level) take into account the 
errors on temperature and /3 model parameters. 
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Fig. 8. Comparison of the different scaled emission measure 
profiles of RXJ0256. 5+0006 with the reference profile (lines) 
based on nearby clusters from Arnaud et al. (2002a). The blue 
(black, red) circles show the NE (overall, SE) profile. The error 
bars on the scaled emission measure are 1 a taking into account 
statistical and temperature uncertainties. 

In the case of hydrostatic equilibrium, the gravita- 
tional mass M to t of a galaxy cluster can be written as : 



M tot (< r) = - 



kT 
Gpm % 



din? 



din 7' 



cllnT 
dlnr 



(8) 



where k is the Boltzmann constant, T the gas tempera- 
ture, G the gravitational constant, p the mean molecular 
weight of the gas {jjl ~ 0.6), m p the proton mass and n e 
the electron density. 

In Sec. 5.1 and 5.2, we found no indication for im- 
portant temperature variations, therefore we neglect the 
temperature gradient in equ.8. For our calculation of the 
mass of the cluster we use in the following temperature 
fcT = 4.9t°jkeV. 

Including the /? model in equ.8 we obtain: 



M tot {< r) 



3k(3 
Gpm p 



T- 



(9) 



Since we have different fit results for the (3 model pa- 
rameters (depending on the region in which the surface 
brightness profile was extracted - see also Sec. 4) we cal- 
culate in the following the mass of the cluster using the 
different model parameters. In order to be able to compare 
our results with other studies we calculate for the different 
(3 models the mass in a region in which the mean density 
of the cluster has an over-density of 500 with respect to 
the critical density of the universe at the cluster redshift 
p c (z). The corresponding radius and mass for over-density 
500 arc defined as rsoo and Afsoo- In a ACDM Universe 
we can write for p c (z): 



Pc(z) = Pco(^m(l + z) 3 + Q A ) 



(10) 



where p c o is the local value of the critical density 
(Weinberg 1972) and : 



PcO 



3flg 

8ttG 



69.4M /kpc 3 



(11) 



At a redshift of 0.36, for fi m =0.3 and Oa=0.7, we find 
p c = lOO.9Af /kpc 3 . 

Using equ.(9) at rsoo we can write : 



M to t(< rsoo) 



M,, 



3fc/3 T 3 
Gjl^; 1 r 500 



r 500 



— 7rr 



.-.in) 



x 500p c 



(12) 



We derive r 50 o using equ.(12) for the different /3 model 
parameters (see Tab. 2). To keep a confidence level of 90%, 
we propagate the corresponding errors from the cluster 
temperature and the (3 model parameters. 

6.4. Comparison of total mass with other work 

We compare our result with the M — T study presented 
by Finoguenov, Reiprich & Bohringer (2001) based on the 
HIFLUGCS galaxy cluster sample (Reiprich & Bohringer 
2002). Adding redshift evolution effects and transposing 
their found relation in our cosmological model we can 
write (neglecting uncertainties): 



-1/2 



M; 



500 



(1 + *) 



-3/2 



k -l/2 
*z=0 



-4.310 13 Mr. 



f kT 



1.58 



(13) 
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-3/2 
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-1/2 
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-0.76 



kT 
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(14) 



Inserting our measured cluster temperature kT = 
4.9 keV, we find M 500 = 3.9x 10 14 M Q and r 500 = 1.2 Mpc. 
These values are in good agreement with the values we 
found with the hydrostatic equilibrium (Tab. 4). Using 
the NE profile, which very likely represents the relaxed 
part of the cluster (see Sec. 6. 2. 2), we found 3.7T^'q x 
10 14 M Q when we use the fit parameters taking into ac- 
count the central part of the surface brightness profile and 
4.1+Qg x 1O 14 M when neglecting the inner parts of the 
surface brightness profile. Discarding evolutionary effects 
(neglecting the (1 + z)~ 3 / 2 term), we find for the above 
relations: r 500 = 2.0 Mpc and M 500 = 6.2 x 10 14 M . 
Not surprisingly, the calculated non-evolution M500 and 
r5oo are much larger than the values we found using the 
best fit (3 model parameters of the NE profile (Tab. 4). The 
value for the no-evolution r§oo exceeds the calculated t^qq 
value even when taking into account the errors. The no- 
evolution M500 agrees within the error bars only with the 
measured M 500 when the central bins are cut out. In this 
case the the calculated uncertainties are substantial (see 
Tab. 4). We take this as an indication that M 500 and rsoo 
do in fact evolve with redshift. 



6.5. Gas mass and gas mass fraction 

We obtain the gas mass by integrating equ. (2) and us- 
ing the /3 model parameters determined in Sec. 4. The 
fraction of X-ray emitting gas of the main cluster of 
RX J0256. 5+0006 is the ratio of gas mass to total mass, 
which we determined in Sec. 6. 3. Our results on the gas 
mass fraction are displayed in Tab. 4. 

The calculated gas masses depend on the determined 
(3 model parameters. Since we find different /3's, not sur- 
prisingly, the different gas masses differ by a factor of 2. 
However, since the errors are relatively large, the gas mass 
fractions agree within the uncertainties. Our found values 
for f g agree within the error bars 

with the value of 20.0+1.9% obtained for hot clus- 
ters from Arnaud & Evrard (1999). For the NE profile 
of the cluster, which seems to represent the relaxed part 
of the cluster (see above and Sec. 6. 2. 2), we find a gas 
mass fraction roughly between 18% and 20%, which is 
slightly lower but in good agreement with the value found 
by Arnaud & Evrard (1999). 

Since the gas mass fraction in galaxy clusters should 
remain constant with redshift, Pen (1997) suggested that 
f g can be used to determine cosmological parameters. If 
we adopt Q m = 1 and J7a = the gas mass fraction lowers 
by about 15%. This results in a gas mass fraction which is 
lower than the value found by Arnaud & Evrard (1999). 
However, since our error bars are large we do not estimate 
this result significant enough to give any statement on the 
value of cosmological parameters. 



Fig. 9. Optical image of the RX J0256. 5+0006 sky region ob- 
tained at the 1.5 m Danish telescope at ESO. The black vertical 
line is a dead column. The contours are the same as in Fig. 2. 
The galaxies with spectroscopic redshift are labeled and each 
of them is tagged by a red circle. The redshifts are presented 
in Tab. 5. 



Label Right ascension Declination 



Redshift 



A 
B 
C 
D 



2^56" 
2 h 56 n 
2 h 56 n 
2 h 56 n 



35.5 s 00° 06' 10.8" 0.36148+0.00051 

33.9 s 00° 06' 40.9" 0.35774+0.00099 

32.8 s 00° 06' 25.3" 0.36451+0.00086 

30.8 s 00° 06' 03.5" 0.37044+0.00076 



Table 5. Spectroscopic redshifts for the labeled galaxies (see 
Fig. 9) obtained at the 4 m Kitt Peak telescope with a 1 a con- 
fidence level (internal errors). 



7. Optical observations of RX J0256. 5+0006 

The SHARC survey collected serendipitous extended 
sources from ROSAT data. An optical follow-up campaign 
of the galaxy clusters was then undertaken (Romer et al. 
2000, Burke et al. 2003). We present in this section the 
optical data obtained for RX J0256. 5+0006. The image 
which was obtained with the 1.5 m Danish telescope at 
ESO is shown in Fig. 9. Redshifts of four galaxies (these 
galaxies present only absorption features, no emission lines 
were present in the optical spectra) were obtained with the 
4m Kitt Peak telescope (see Tab. 5 and Fig. 9). 

The error bars obtained for the galaxy redshifts in 
Tab. 5 only account for the internal errors due to cross 
correlation analysis. The redshifts of galaxies B and D ob- 
tained from the SDSS (Sloan Digital Sky Survey) are 0.358 
and 0.371, respectively. These values are in good agree- 
ment with the observed redshifts in Tab. 5 and indicate 
that the SHARC redshifts arc correct. 



7.1. Galaxy position and X-ray emission correlations 

Fig. 9 displays many bright galaxies inside the region of 
the main X-ray cluster. 

One galaxy, located approximately 8" south of galaxy 
B coincides nicely with the central residual found in X- 
rays (see Sec. 3). Since the extent of this residual is very 
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small, it is not unlikely that the corresponding source is 
point-like, such as an AGN for example. 

Galaxy D coincides well with the X-ray maximum of 
the subcluster. This galaxy is by far the brightest in this 
region. There is an offset of roughly 5" between galaxy D 
and the maximum of the residuals. The projected offset 
corresponds to 40kpc. An explanation for this offset can 
be ram pressure stripping, which acts on the ICM but not 
on the galaxies during infall: unlike galaxies, which can 
be considered as collisionless matter particles which follow 
directly the dark matter distribution in cluster mergers, 
the hot gas of the infalling subcluster is slowed down due 
to presence of the ICM of the main cluster (see numerical 
simulations performed by e.g. Roettiger, Loken & Burns 
1997; Stevens Acreman & Ponman 1999). This effect was 
for example studied in Neumann et al. (2001) in the case 
of the Coma cluster. We will address the issue of ram 
pressure stripping in detail in Sec. 8. 3. 

7.2. Velocity differences between main and subcluster 

In the following we will estimate the velocity difference 
along the line-of-sight v± between the main cluster and the 
subcluster component based on the redshift measurements 
of the four galaxies (see also Tab. 5). xSince we only have 
the redshift for one galaxy in the subcluster, galaxy D 
(z = 0.3704) we assume that this redshift is representative 
for the entire subcluster. 

Calculating the mean redshift of the other three galax- 
ies (A, B, and C), assuming that they all belong to the 
main cluster and neglecting peculiar motions, we find 
z = 0.3612. The difference in redshift between this mean 
redshift and the redshift of galaxy D corresponds to a ve- 
locity of 2750 km/s. In order to estimate the minimal v± 
we compare the redshift of galaxy D with the redshift of 
galaxy C (z — 0.3645), which has of the three galaxies A, 
B and C the redshift closest to galaxy D. The velocity dif- 
ference corresponds in this case to 1750 km/s and is thus 
the minimal v±. 

Since there exists very likely a velocity component par- 
allel to the plane-of-sight, v± is only a lower limit on the 
actual infall velocity. In the following we will estimate the 
infall velocity and calculate a possible merger model for 
this cluster. 



8. Dynamics of the merger in RX J0256. 5+0006 

8.1. Geometry of the merging cluster components 

In the following we will try to quantify the merger geome- 
try in this cluster by using the estimates of the projected 
distance between the two clusters d m i n (350 kpc — see also 
Sec. 3) and vj_ (Sec. 7. 2). The adopted merger geometry is 
shown in Fig. 10. 

A similar approach has been already performed by e.g. 
Beers, Geller & Huchra (1982) for Abell 98 and Colless & 
Dunn (1996) for the Coma cluster. 




Earth 

Fig. 10. Merger geometry of RX J0256. 5+0006 in the case of 
an on-axis merger. The small point in the centre represents the 
main cluster centre, and the larger point represents the sub- 
cluster. d m in is the projected distance between the two clus- 
ter centers (see Sec. 3), d is the physical distance between the 
two components) and v± , the projection of the impact velocity 
along the line-of-sight (see Sec. 7. 2). We estimate that the main 
cluster extends up to r2oo- 9 is the impact angle. 



8.1.1. The merger model 

To determine the merger geometry analytically, we apply 
the following assumptions : 

- a purely gravitational two-body interaction, i.e. pressure 

free and no neighborhood effects ; 

- the main component is defined by the mass distribution 

M t ot{d) (see equ. (9)) ; 

- the subcluster mass distribution can be approximated 

with a point ; 

- The main cluster is at rest, only the subcluster moves. 

From the law of gravitation follows : 
dv _ GM tot {d) . 



dt 



d 2 



(15) 



v is the velocity of the subcluster, d is the distance be- 
tween main cluster center and subcluster. i is an unity 
vector along the merger axis and is oriented in the direc- 
tion of increasing d. 
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Fig. 11. The velocity of the infalling subcluster as function of 
its distance d to the main cluster centre. The full line represents 
the velocity when using the mass profile of the /3 model param- 
eters from the NE surface brightness model (see also Tab. 2). 
The dotted line shows the velocity when using for the mass 
distribution of the main cluster the /3 model parameters from 
the overall surface brightness profile with the subcluster cut 
out. 

Using for simplicity that the distance of the two cluster 
components was close to infinity at the beginning we can 
distinguish between two scenario : 

- d > r 2 oo : (numerical simulations suggest that r2oo corre- 
sponds to the virialized part of a cluster) in this case, 
the main cluster can be approximated with a point 
mass located at the cluster center. The velocity is : 



2GM 2 



(16) 



d < f2oo : in this case, we have to take into account the 
mass distribution defined in equ. (9). By integrating 
equ. (15), we can write : 



\ 



2GM 20 o 3fc/?T / d 2 +r 2 



7*200 



/imp 



(17) 



'200 



To determine r2oo we apply equ. (12) with an over- 
density of 200. For the main cluster mass distribution as 
well as r2oo and M200, we use the results of the (3 model 
fits. More specifically we use in the following the fit results 
of the NE surface brightness profile and the overall profile 
taking into account the central part of the cluster (see 
Tab. 2 - the a case). The comparison of the results based 
on these two models gives an idea on the dependence of 
the merger dynamics on the mass profile. We give more 
weight on the results of the NE profile, since it seems to 
represent the relaxed part of the main cluster in which the 
hypothesis of hydrostatic equilibrium is most likely valid. 

Using the overall model without angular selection we 
find M 200 = 9.2 x 10 14 M Q and r 200 = 2.2 Mpc. For the NE 
model we find M 20 o = 5.9xlO 14 M and r 20 o = 1.9 Mpc. 

Fig. 11 shows the obtained relationship between dis- 
tance d and velocity v for the two models. 
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Fig. 12. The vertical velocity component v± against the clus- 
ter distance according to the model described in Sec. 8. 1.1 and 
equ. (20). Full and dotted line follow the same convention as 
Fig. 11 



8.1.2. Impact angle and physical distance 



From Fig. 10 follows 



i9 = 



d„ 



cosf = 



d 

v± 



(18) 
(19) 



9 is the impact angle, which is defined as the angle between 
the line-of-sight and the merger axis. The merger axis is 
the line which connects the two cluster centers. 

Since we observe a comet-like structure for the sub- 
cluster (which indicates gas stripping, see also below), we 
assume that the two cluster components are interacting 
which implies d < r2oo- 

The redshift of galaxy D (see Tab. 5) which we assume 
to be located roughly at the center of the subcluster is 
greater than the redshifts obtained inside the main cluster. 
This implies that the subcluster is in front of the main 
component and thus that 9 must be lower than 90° (see 
Fig.10). 

The impact geometry can be furthermore restricted by 
combining equ. (18) and equ. (19) to: 



i 2 ■ 

nun 

d 2 



% = 1 

V 2 



(20) 



The relationship between v and d is shown in Fig.ll. 
Since we do not possess a precise estimate of v±, we 
show in Fig. 12 the relation between v± and d. We can 
set the following limits: d m i n < d < r2oo (Sec. 3) and 
1750 km/sec< v± <2750 km/sec (Sec. 7. 2). To constrain 
v± more precisely we would need more galaxy redshift 
measurements. 

Fig. 13 and Fig. 14 show the dependence of the impact 
angle 9 on cluster distance d and v±. The greater 9 is, 
the smaller is d. Based on our constraints on vj_ (vertical 
line in Fig. 14), we obtain values for 9 between 15°and 
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Fig. 13. The impact angle 6 defined in Fig. 10 as function of 
distance d between the two merger components. Since the pro- 
jected distance is 350 kpc, d > 350 kpc. - Please note: the 
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Fig. 14. Impact angle 9 defined in Fig. 10 versus v±. Full and 
dotted line follow the same convention as Fig. 12. The vertical 
line shows the minimal value for v±. 



44°(9°and 60°) for the NE (overall) mass model. Fig. 14 
shows that for each of the two adopted mass models any 
given value of v±, gives two possible solutions for 9. This 
means that a measured value of v± gives two different 
values of 9 connected to two possible physical distances d. 
For 15° < 6 < 44° wc find for the distance d: 0.5 Mpc< 
d <1.4 Mpc for the NE model. 

8.1.3. Merger model dependence on main cluster mass 
distribution 

We see that the smaller the total mass of the main cluster 
is, the tighter are the constraints on the geometry of our 
merger model and angle 9. This is not surprising since the 
total infall velocity is a function of cluster mass. Smaller 
masses show generally smaller infall velocities. For the NE 
mass model the observed vertical velocity v± is close to the 



maximal possible vj_ of about 1900 km/s. For comparison, 
for the overall mass model case the maximum velocity is 
roughly 2500 km/sec. 

8.2. The cometary shape of the subcluster: effects of 
ram pressure stripping 

The subcluster has a clear cometary shape, which suggests 
that ram pressure due to the infall plays an important role. 
In the following we will estimate the importance of ram 
pressure stripping in the merger model. Ram pressure can 
be denned in our case as: 



Prp — 



(21) 



p g is the gas density of the main cluster (p g = 
np,m p ; \i = 0.6; m p : proton mass) and v is the infall ve- 
locity of the subcluster with respect to the main cluster 
centre. The calculation of the merger geometry suggests 
that the subcluster has a physical distance of the main 
cluster of 0.5 Mpc< d <1.4 Mpc. Effects of ram pressure 
are important if it is of the same order or larger than the 
internal pressure of the subcluster. The internal pressure 



of the subcluster p 



«<> 



ikT = 7 x 10~ 12 g cnT 



suming an electron number density of n e (subcluster) = 
7.1 x 10 -4 cm -3 (see Sec. 8. 5) and an internal temperature 
of the subcluster of about 3 keV (3 keV is the estimated 
temperature of the subcluster when using the z-evolving 
Lx — T-relation and the measured Lx = 1.9 x 10 44 erg/s 
of the subcluster found below). Using as a typical infall 
velocity of v = 2000km/s we find that p g must be larger 
than 1.8 x 10~ 28 g cm~ 3 so that p rp > p sg . Using the NE 
(overall) /3-modcl parameters wc find that for example at 
1 Mpc: p g = 4.5 x 10~ 28 gcm~ 3 (p g = 5.5 x 10~ 28 gcm~ 3 ). 
Since these values are larger than 1.8 x 10~ 28 gcm~ 3 we 
conclude that ram pressure is actually important at this 
merger state. We display in Fig. 15 the ram pressure as 
function of distance to the main cluster centre for the over- 
all and the NE profile. We use equ.(17) as velocity input 
for the calculation of the ram pressure in equ.17 (see also 
Fig. 11). We equally show the value of the internal pressure 
of the subcluster assuming a constant gas density distri- 
bution as horizontal line in Fig. 15 for comparison. 

In our adopted merger geometry, in which we find for 
the impact angle 15° < 6 < 44°, the physical length of 
the subcluster is much larger than its projected length. 
If the subcluster for example had the form of a line then 
for 9 — 20°, which corresponds to a physical distance of 
1 Mpc the real length of the subcluster was three times 
larger than its projected length. The observed length is 
approximately 700 kpc, which implies a physical length 
in case of a line of approximately 2 Mpc (this is a factor 
of 3 larger than the width -600 kpc- of the structure). 
However, it is very unlikely that the subcluster has the 
form of a line, it is more likely for it to be ellipsoidal. In 
this case the calculation of the physical length of the sub- 
cluster is more complicated. In the following we attempt 
to calculate the length of the subcluster assuming an el- 
lipsoidal shape. Fig. 16 shows our adopted model for the 
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Fig. 15. The ram pressure as a function of distance to the main 
cluster centre using equ.17 and the gas density distribution 
based on the j3 model parameters. The full line corresponds to 
the NE and the dotted line to the overall f3 model parameters. 
The horizontal dashed line indicates the internal pressure of 
the subcluster if we assume a constant electron gas density of 
n e - 7.1 x 10~ 3 cm~ 3 (see also Sec.8.5). 




Fig. 16. The ellipse model of the subcluster from above, look- 
ing down from North to South. 2x is the observed projected 
length of the subcluster. 



ellipsoid from above looking downwards (from North to 
South). 

We use for the calculation of the extent of the subclus- 
ter the ellipse equation: 



b 2 



1 



(22) 



in which x' and y' are the x and y coordinates before 
rotation around the impact angle 9. We furthermore use 
for the rotation around 9: 



x' sin 9 + y' cos 9 = x sin 9 



Where x is the x-coordinate after rotation, ext (see 
Fig. 16 is the maximum value in the x-direction after ro- 
tation. We can therefore write: 



d(x = ext) 
dx' 







Using equ.23 and equ.24 we find: 




h 2 



i/2 V " V * "ea^tan 2 ^ 
Inserting equ.25 in equ.23 we find 



ext = x sin ( 



\ 



1 - 



1 



1 



Ab 2 
x' 2 tan 2 I 



(24) 



(25) 



(26) 



1 -6cos( 

a z 



(23) 



Unfortunately this equation is too complicated to cal- 
culate x' directly as function of 9. However, knowing ext 
and b (we assume symmetry around the merger axis, so 
b = 300 kpc and ext = 350 kpc see also Sec. 3) we can 
search for a given of 9 the corresponding x' with the re- 
quirement that exi=350 kpc. Subsequently we can calcu- 
late a(x'). 2a is the physical extent of the subcluster along 
the merger axis, c = 2a + d is the distance of the outer 
boundary of the subcluster to main cluster centre (d and 
8 are correlated according to Fig. 13). From Fig. 15 we can 
give limits on c. Applying the internal pressure argument 
of the subcluster from above we find a maximum value for 
c of 1.5 Mpc. If we assume that r2oo defines the area in 
which ram pressure plays an important role then we find 
a maximum of c = 1.9 Mpc. Fig. 17 shows c as function 
of distance to the cluster centre. We find that distances 
up to 0.9 Mpc are possible in case of r2oo and 0.6 Mpc in 
the case of c = 1.5 Mpc. If we assume for example that 
2b = 400 kpc instead of 600 kpc, we find the results indi- 
cated as dotted line in Fig. 17. In this case we find smaller 
values for d. This shows the dependence of our calcula- 
tions on the choice of b. The smaller b, the smaller is the 
allowed range of d. However, it seems rather unlikely that 
2b = 400 kpc, which would imply a relatively large asym- 
metry of the subcluster. Ram pressure stripping makes 
the substructure more symmetric around the axis of in- 
fall (the parts of the subcluster with the largest distance 
with respect to its centre are less bound to the structure 
and are thus more easily stripped). Since ram pressure is 
efficient in this merger it is very likely that the symmetry 
around the infall axis is already established. 

We conclude that due to observed effects of ram pres- 
sure stripping that the subcluster has a distance to the 
main cluster of roughly 0.6-0.9 Mpc. 

8.3. The offset between X-ray maximum and galaxy D 

We have seen that ram pressure stripping is efficient in 
the subcluster as it falls onto the main cluster. Ram pres- 
sure acts mostly on the gas and much less on galaxies. 
Therefore there should be a displacement between the 
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Fig. 17. The subcluster outer boundary with respect to the 
main cluster centre (c = 2a + d) as function of subcluster main 
cluster distance d (see also text). 

subcluster ICM and galaxy D, which we assume traces 
the centre of gravity of the subcluster (which is supposs- 
edly dominated by dark matter). We actually see a dis- 
placement of about 5" (see Fig. 2). We can estimate the 
expected displacement between gas and galaxy D using 



Mg as vdec = / M gas decdt = 



PrpS—dr = 



Prp 



dr(27) 



with Mgas'- gas mass of subcluster; Vdec the difference of 
velocity between galaxy D and gas; dec: deceleration of 
the gas with respect to the subcluster centre; S: Surface 
of subgroup facing the centre of the main cluster. We as- 
sume that the gas density is constant within the subclus- 
ter. We integrate equ.27 numerically. We use for M gas 
10 13 M Q , which is about 10% of the total mass found for 
the subcluster (see Sec. 8. 5). For S we use a circle with ra- 
dius 300 kpc (see above), which likely matches half of the 
diameter of the structure in width. Vd ec is inversely pro- 
portional to Mgas and directly proportional to S. In order 
to estimate the displacement of the centre of the gas with 
respect to galaxy D, we integrate 



Vdecdt = 



dt , 
Vdec-rdr 
dr 



I'de 



-dr 



(28) 



Again, we integrate numerically to calculate ddi sp the 
distance between galaxy D and the location of the X-ray 
residual. In Fig. 18 we show ddisp as a function of the dis- 
tance of the subcluster to the main cluster centre. Fig. 19 
shows the projected distance of galaxy D and subclus- 
ter in arcmin as function of main cluster distance. The 
full (dotted) line corresponds to the NE (overall) mass 
model. The horizontal dashed line corresponds to the dis- 
placement observed in Fig. 9 of approximately 5" , which 
implies a distance of subcluster to main cluster centre of 
about 700 kpc. It is difficult to assess the correct offset 
between the X-ray maximum of the subcluster and galaxy 
D. Several parameters can in principle play a role: the 
parameters of the subtracted (3 model and uncertainties 



Fig. 18. The physical distance between galaxy D and the max- 
imum of X-ray emission of the subcluster based on calculations 
of ram pressure stripping as function of distance to the main 
cluster centre. Full line: calculations based on NE /3 model ; 
dotted line: based on overall /3 model . 



of the alignment of X-ray image and optical image. We 
vary the parameters of the subtracted model within the 
error bars and see that the location of the X-ray residual 
maximum does not change. We estimate the alignment un- 
certainties, which are correlated with PSF effects in the 
order of half of the PSF size of Xmm-NEWTON , which is 
about 2-3 arcsec. Therefore, we estimate the maximum al- 
lowed offset between X-ray maximum and galaxy D to be 
8 arcsec, or 0.133 arcmin. The implied minimal distance 
of the subcluster to the main cluster centre in our simple 
model is thus about 600 kpc (see also Fig. 19). Inversely, 
the minimal offset between galaxy D and the subcluster 
X-ray emission is about 2 arcsec. Using Fig. 19, this corre- 
sponds to a distance of the subcluster to the main cluster 
of about 1 Mpc. 

Therefore, the distance of the subcluster to the main 
cluster calculated from a simple ram pressure stripping 
model varies between 0.6 to 1 Mpc. This suggests that 
the impact angle 8 lies between 20 to 35°. This is in good 
agreement with our merger model. We would like to stress, 
however, that the model we adopted for ram pressure is 
very simple and is based on many assumptions. 

8.4. Constraints on impact parameter 

We find coherent results for the distance of the subcluster 
to the main cluster centre using our on-axis merger model 
and calculating effects of ram pressure. With all this infor- 
mation we are in principle able to put some constraints on 
the probability that the merger is off-axis with a non-zero 
impact parameter. Most possible scenarios for an off-axis 
merger would increase the visible effects of ram pressure 
stripping since off-axis merger trajectories are generally 
longer than on-axis merger trajectories. Therefore there is 
more time for ram pressure to be efficient. We would thus 
expect that ram pressure calculations give distances of the 
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Table 6. Corrected luminosity L x and mass M2 of the subclus- 
ter depending on the distance of the subcluster to main cluster. 
The main cluster density n e ,i is calculated according to the /3 
model presented in Sec. 4. The mass and the temperature T 
of the subcluster are estimated according to the scaling laws 
presented in Sec. 6. 3. The mass ratio between main and sub- 
cluster is computed assuming for the main cluster component 
M200 = 5.9 x 10 14 M Q (see also Sec.8.1.1). 



Fig. 19. The projected distance between galaxy D and the 
maximum of X-ray emission of the subcluster based on cal- 
culations of ram pressure stripping as function of distance to 
the main cluster centre. Full line: based on NE /3 model ; dot- 
ted line: based on overall /3 model . The horizontal dashed 
line shows the observed offset of galaxy D to subcluster X-ray 



maximum. 



subcluster to the main cluster, which are smaller than dis- 
tances calculated with the on-axis merger model. However, 
this is not the case here. Remaining possible models are 
those in which beside the infall velocity related to the grav- 
itational force there exists an additional velocity compo- 
nent in direction of the line-of-sight and perpendicular to 
it with a specific relationship between the velocity com- 
ponents. In this particular case, by chance, the merger 
would have a non-zero impact parameter with a specific 
orientation, which mimics at the moment of observation 
an on-axis merger with the correct projected elongation 
of the substructure caused by ram pressure. It is doubtful 
whether this is likely. A simple explanation of an on-axis 
merger seems to us much more realistic. 

8.5. Subcluster mass estimate 

Determining the luminosity of the subcluster is more diffi- 
cult than for the main cluster (see Sec. 6.1) since we cannot 
fit an adequate model to the comet-like structure. To over- 
come this problem, we simply determine the count rate in 
the region of the residual map (after main cluster emis- 
sion subtraction, see also Sec. 3) in which the subcluster is 
clearly detected. Translating the count rate into luminos- 
ity, we find Lboi, x = 1.9 ± 0.3 x 10 44 ergs/s. 

Since we use the luminosity estimate via the scaling 
law relations to obtain the mass of the subcluster, lumi- 
nosity biases introduce biases in our mass estimate. The 
luminosity estimate of the subcluster can be biased so 
that : 

- we underestimate the total count rate since we only take 
into account the region in which we detected the sub- 
cluster. This subcluster might be extended to larger 
regions for which we do not account. 



since the bremsstrahlung emission evolves with r? e and 
because the subcluster lies very likely inside the main 
cluster gas sphere, the main cluster emission model 
subtraction is not correct to obtain the subcluster 
emission. In the case that the subcluster is in the gas 
sphere of the main cluster, the total emission evolves 



with 



(n eA +n e , 2 y >n 2 el +i 



where 



te,l 



is the 



electron density of the main cluster, and n e ^ the elec- 
tron density of the subcluster. The quantity 2n e _in e _2 
enhances the measured luminosity. The value we mea- 
sure by subtracting the main cluster component is the 
luminosity of the subcluster in the case where it is 
not physically connected to the main cluster. We thus 
might overestimate the luminosity of the subcluster. 

We concentrate in the following on the second case, i.e. 
correct the estimated luminosity for the remaining main 
cluster contribution. 

We assume axis symmetry around the merger axis of 
the subcluster. The minor axis — about 1.4' or 600 kpc 
, measured in the N-S direction (see also Sec. 3) is as- 
sumed to be equal to the minor axis of the subcluster 
perpendicular to the N-S axis and to the merger axis. 
Furthermore we assume that the subcluster electron den- 
sity is constant over its volume. If we furthermore assume 
d = d m i n = 350 kpc we can calculate n e _i by integrating 
the )3 model representing the main cluster density distri- 
bution (Sec. 4) along the line-of-sight. In this case we find, 
applying rig = (n ey i+n ey 2) 2 '- H e ,2 =7.1x 10~ 4 cm -3 for the 
subcluster electron density. 

Since we cannot determine the physical distance d of 
the subcluster to the main component with high accuracy, 
we show in Tab. 6 the luminosity and the mass of the sub- 
cluster and the mass ratio between the two components 
for several values of d (0.5 < d < 2 Mpc). 

Since the subcluster has very likely a physical distance 
in the order of 0.5-1 Mpc (Sec. 8. 1.2), we estimate the mass 
ratio M 8U b/M ma i n to be between 20-30%. This implies 
that RX J0256. 5+0006 is a major merger event. 
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8.6. The observed subcluster temperature and 
implications 

One important question is why the region of the subcluster 
shows such a high temperature. According to the luminos- 
ity estimated in Sec. 8. 5 and equ. (5), the highest temper- 
ature expected for the subcluster would be about 3.2 keV. 
However, we measure a temperature of 5.7 _{ keV in this 
region. 

It is clear that the emission of the main cluster con- 
taminates the subcluster region. However, if this was the 
only contamination, the measured temperature should lie 
between 3 to 5 keV, but not above 5 keV. This suggests 
the presence of a third temperature component in this re- 
gion. This component is quite likely the region in front of 
the subcluster pointing towards the main cluster centre, 
which is heated cither by adiabatic compression or shock 
waves due to the infall. The merger geometry, which favors 
small values for the impact angle 9, maximizes the effect 
of projection, since the heated region and the subcluster 
region lie in this case roughly in the same line-of-sight. 

In order to estimate the temperature of this third hot 
component we perform a three temperature fit to the spec- 
trum of the subcluster region. We fix the temperature of 
the main cluster component (determined from our spec- 
tral fits) as well as the subcluster component (determined 
by converting the observed luminosity of this component 
in temperature). The intensity of the main cluster compo- 
nent is fixed and calculated from the j3 model parameters 
and the subcluster emission is fixed and estimated from 
the count rate estimate after main cluster contribution. 
Unfortunately the statistics of the spectrum is too poor 
to allow the temperature determination of this third hot 
component. 

9. Discussion and conclusion 

We present results of an Xmm-NEWTON observation of 
RX J0256. 5+0006 a medium distant cluster of galaxies 
(z=0.36) found in the Bright SHARC catalog. We observe 
two X-ray maxima in this cluster. The principal X-ray 
maximum is linked to the main cluster centre and a sec- 
ond X-ray maximum lies West with respect to the main 
cluster centre. The observed bimodal morphology of the 
cluster indicates that this cluster is in a merger state. 
Hydrodynamic simulations have shown that merging ac- 
tivity in clusters is naturally linked to temperature gradi- 
ents in the ICM. In our spectro-imaging analysis, however, 
we do not find strong indications for temperature varia- 
tions. The lack of temperature gradients is confirmed by 
spectral fitting in selected regions. 

For the main cluster centre we find a temperature 



of kT 



4.9 



-0.6 
-0.5 



keV, which is representative for the 



entire main cluster component. We determine the bolo- 
metric X-ray luminosity of the main cluster to be Lx = 
1.5 x 10 45 erg/s. This luminosity is higher than the value 
(Lx ~ 9x 10 44 erg/s) expected from the luminosity- 
temperature relation established by Arnaud & Evrard 



(1999) taking into account redshift evolution (the term: 
(1 + z) 3 / 2 ). Neglecting redshift evolution the discrepancy 
between observed and expected luminosity becomes even 
stronger. We take this as an indication that X-ray lumi- 
nosity evolves with redshift, as predicted from simple the- 
ory on structure formation. Our hint for evolution agrees 
with a study by Vikhlinin ct al. (2002) based on clusters 
with (z > 0.4). The high X-ray luminosity of this cluster 
is probably also linked to the fact that the main cluster 
in itself is not entirely relaxed. This latter hypothesis is 
supported by the fact that the surface brightness profiles 
extracted in different sectors show considerable variation 
in shape. 

X-ray profiles of clusters show a remarkable degree of 
self-similarity (see for example Neumann & Arnaud 1999). 
We compare the X-ray profiles of RX J0256. 5+0006 with 
the reference profile based on nearby clusters defined in 
Arnaud, Aghanim & Neumann (2002). We see that only 
the surface brightness profile extracted North-East (NE) 
of the cluster agrees with the reference profile of Arnaud 
et al. (1999). The other profiles (the overall profile with 
the subcluster cut out and a profile extracted SE of the 
cluster) show larger core radii with respect to the refer- 
ence profile. We take this as indication that the NE sector 
represents the relaxed part of the cluster. 

Fitting a j3 model to the different surface brightness 
profiles, we observe that r c (between 100 and 550 kpc) 
and (0.57-1.2) vary strongly between the different pro- 
files. Consequently the mass estimates based on hydro- 
static equilibrium and (3 model parameters are quite dif- 
ferent. We observe M 500 between 3.7 and 12 x 10 14 M Q . 
For the NE profile, which seems to represent the relaxed 
part of the main cluster we find M500 ~4x 10 14 M© . 
This is in good agreement with a z-evolving M — T re- 
lation, when compared to the results of the HIFLUGCS 
sample (Finoguenov, Rciprich & Bohringcr 2001) which 
gives M500 = 3.9 x 1O 14 M0 (transposed into our cos- 
mological model) . A non-z-evolving M — T relation gives 
-M500 = 6.2 x 1O 14 M0 , which is only marginally consistent 
with our observed mass for the main cluster based on the 
NE surface brightness profile. 

With the f3 model parameters we also calculate the 
gas mass of the cluster. The corresponding gas mass frac- 
tion (f g ) of the cluster varies between 12 and 33%. Taking 
into account only the overall and the NE profile we find 
fg's in the order of 18-20%, which is in good agreement 
with the value found by Arnaud & Evrard (1999). Pen 
(1997) suggested that constant f g in clusters can be used 
to determine cosmological parameters. Changing cosmol- 
ogy to Q m = 1, we find indeed lower f g 's around 16- 
18%. This is a small indication that indeed Vt m = 0.3 and 
Oa = 0.7, as used in this paper. However, Neumann & 
Arnaud (2001) found that the gas mass in clusters does 
not follow M gas ex T 3 / 2 , but M gas ex T 194 . If M tot ex T 3 / 2 
as predicted from simple structure formation, or generally 
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M to t oc T x with x 7^ 1.94 2 , then f g varies as function 
of cluster temperature and indicates that internal cluster 
physics plays an important role. In this case, using gas 
mass fractions in clusters for the determination of cosmo- 
logical parameters is more complicated. One needs to com- 
pare f g of nearby and distant clusters at exactly the same 
temperature. This approach is furthermore only valid if 
additionally there does not exist an intrinsic redshift de- 
pendence of f g (see for example Ettori et al. 2004 for the 
redshift dependence of several physical quantities in clus- 
ters). If internal cluster physics change significantly f g in 
clusters it is very likely that this effect is redshift depen- 
dent. If this is true than using f g for determining cosmo- 
logical parameters depends on a large number of physical 
relationships, which would make this approach less attrac- 
tive. 

Beside the properties of the main cluster of 
RX J0256. 5+0006 we also study the cluster's merger 
dynamics. We develop a simple model to constrain the 
merger geometry based on an on-axis merger with impact 
parameter zero. As input we use the projected distance 
of the subcluster to the main cluster centre observed in 
X-rays (350 kpc) and estimates of 4 cluster galaxy red- 
shifts obtained with the Kitt Peak telescope. The possible 
velocity difference between main and subcluster lies be- 
tween roughly 1800 and 2800 km/s. The range of possible 
distances ranges between 0.5 Mpc< d < 1.4 Mpc, which 
limits the impact angle 6 to: 15° < 6 < 44°. The rela- 
tively large distances of the subcluster to the main cluster 
centre are confirmed by the fact that we do not see tem- 
perature variations along the merger axis (see above). We 
also examine the effects of ram pressure stripping in the 
subcluster ICM. We find that the effects we see in X-rays 
like cometary shape of the subcluster ICM as well as dis- 
placement of X-ray maximum with respect to the main 
galaxy of the subcluster (galaxy D) suggest distances to 
the main cluster centre of about 0.6 to 1 Mpc, which im- 
ply 20° < 8 < 35°. This is in very good agreement with 
the simple merger model and suggests that the impact 
parameter of the merger is very likely close to zero. 

Estimating the mass of the subcluster using a 
luminosity- mass relationship we find that the ratio of sub- 
cluster to main cluster mass is in the order of 20 to 30%. 
This ratio is relatively high when compared to other clus- 
ter mergers and indicates that we observe a major merger 
in RX J0256. 5+0006. 
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